Multipartite Entanglement Signature of Quantum Phase Transitions 
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We derive a general relation between the non-analyticities of the ground state energy and those 
of a subclass of the multipartite generalized global entanglement (GGE) measure defined by T. R. 
de Oliveira et al. [Phys. Rev. A 73, 010305(R) (2006)] for many-particle systems. We show that 
GGE signals both a critical point location and the order of a quantum phase transition (QPT). We 
also show that GGE allows us to study the relation between multipartite entanglement and QPTs, 
suggesting that multipartite but not bipartite entanglement is favored at the critical point. Finally, 
using GGE we were able, at a second order QPT, to define a diverging entanglement length (EL) 
in terms of the usual correlation length. We exemplify this with the XY spin-1/2 chain and show 
that the EL is half the correlation length. 
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Quantum Phase Transitions (QPTs) occur at zero tem- 
perature and are characterized by non-analytical changes 
in the physical properties of the ground state of a many- 
body system governed by the variation of a parameter 
A of the system's Hamiltonian H{X). These changes are 
driven solely by quantum fluctuations and are usually 
characterized by the appearance of a non-zero order pa- 
rameter Since QPTs occur at T = 0, the emerging 
correlations have a purely quantum origin. Therefore, it 
is reasonable to conjecture that entanglement is a cru- 
cial ingredient for the occurrence of QPTs (e.g. Refs. 
0) B S and references therein). If this is true, then 
the QPT would imprint its signature on the behavior of 
an entanglement measure. Under a set of reasonable gen- 
eral assumptions, Wu et al. have demonstrated that a 
discontinuity in a bipartite entanglement measure (con- 
currence [l3| and negativity [lH) is a necessary and suf- 
ficient indicator of a first order quantum phase transition 
(IQPT), the latter being characterized by a discontinuity 
in the first derivative of the ground state energy. Further- 
more, they have shown that a discontinuity or a diver- 
gence in the first derivative of the same measure (assum- 
ing it is continuous) is a necessary and sufficient indicator 
of a second order QPT (2QPT), which is characterized 
by a discontinuity or a divergence of the second deriva- 
tive of the ground state energy. Nevertheless, most of the 
models of 2QPTs considered so far did not present any 
long-range bipartite entanglement at the critical point, 
even though the correlation length diverges. Moreover, 
contrary to expectations, most of the measures discussed 
in the literature are, to the best of our knowledge, not 
maximal at the critical point, the exceptions being the 
one-site von Neumann entropy of the Ising chain [3] , the 
localizable entanglement of a finite Ising chain with 14 
sites [l2|, and some classes of the Generalized Global 
Entanglement (GGE) for the Ising chain Q. 

In this Letter we firstly extend Wu et al. M results 
to a multipartite entanglement (ME) measure [8,0], the 
GGE introduced in Refs. [1, Uj], and discuss how non- 
analyticities in the energy are signaled by the GGE. Sec- 



ondly, we define an entanglement length (EL) for an ar- 
bitrary collection of two-level systems. In the case of a 
symmetry-breaking 2QPT, this EL diverges at the criti- 
cal point and is simply related to the correlation length. 
This result indicates that ME is most favored at that 
point, contrary to bipartite entanglement 0, We 
consider the consequences of this result for specific spin- 
1/2 models presenting IQPT or 2QPT. 

In particular, for the 2QPT of the one-dimensional 
transverse field XY model we obtain all the relevant criti- 
cal exponents, with the EL defined in terms of correlation 
functions (CFs) appearing in the GGE. We also show in 
this specific case that the GGE is maximal at the crit- 
ical point, thus signaling the QPT,_as three of us had 
already observed in the Ising case [a]. This last result, 
together with a diverging ME length at the critical point, 
reinforces that ME plays a significant role in QPTs. 

Following Ref. , a discontinuity in (discontinuity in 
or divergence of the first derivative of) the concurrence 
or negativity is both necessary and sufficient to signal 
a IQPT (2QPT) for systems of distinguishable particles 
governed by up to two-body Hamiltonians. The energy 
per particle (e) derivatives depend on the two-particle 
density matrix elements as 

d^e = (l/iV)^TrPAC/(z,j))py], (1) 
die = {l/N)Y,{Tr[{^lU{^,J))p.,] 

+ Tr[(d,.U{i,j))d^p,,]}, (2) 

where pij is the reduced two-particle density operator 
and U {i,j) includes all the single and two-body terms of 
the Hamiltonian associated with particles i and j. Now, 
assuming that U (i, j) is a smooth function of the Hamil- 
tonian parameters and that pij is finite at the critical 
point, the origin of the discontinuity in the energy (dis- 
continuity in or divergence of the first derivative of the 
energy) is the fact that one or more of the elements of 
Pij (d\Pij) are discontinuous (divergent) at the transi- 
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tion point A = Ac Since the concurrence and the 
negativity are both hnear functions of the elements of 
Pij it turns out that a discontinuity/divergence in one 
of them (in the derivative of one of them) impHes a dis- 
continuity/divergence of the energy (in the derivative of 
the energy) and vice- versa 0]- A natural question then 
arises: Does a ME measure show the same feature? In 
what follows, we give an explicit affirmative answer to 
this question 

In [a] three of us introduced two new quantities, both 
of which can be seen as generalizations of the Meyer- 
Wallach lj| Global Entanglement, originally defined 
for a system of N parties (particles). The first one 
is the average linear entropy of all Ni < N parti- 
cles, where we assume a fixed "distance" between the 
A'^i particles. The second quantity is an average over 
all possible distances/configurations in which the A^i 
particles can be arranged [a, For iVi = 1 both 

quantities are the same (G(l) — -Eq"*) and we re- 
cover the Meyer- Wallach measure. The first non-trivial 
case appears when TV > 4 and we pick two parti- 
cles (iVi = 2) labeled by i and j. Now for a den- 
sity matrix Pjj+n of dimension d we have G (2, n) = 



d 
d-l 



1 



which is the 



mean linear entropy of all pairs of particles n = \i — j\ 
sites apart, i.e., the mean entanglement between these 
pairs and the remaining N — 2 particles. Averag- 
ing over all possible distances 1 < n < N, 



(2) 



N(N-i) Yln=i(^ ^ (2: '^)- In order to simplify the 
notation (with no loss of generality), from now on we will 
work with the linear entropy of a single pair of particles 
n sites apart, which we call Q{2,n). Note that in this 
notation G(2,n) = Q{2,n), being the average of 0{2,n) 
over all particles n sites apart. For a translationally sym- 
metric system G(2,n) = Q{2,n). 

Considering Q (2, n) as a function of the tuning param- 
eter A we can write it and its derivative in terms of the 
Im elements of pij ([Pj j+n];„) as 



e(2,n) 



d-l 



1" H 



l.m—l 



(3) 
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dxgi2,n) = -—^J2\[pj,j+n]iJdx\[pjj+„],J.{i) 

Lrn — l 

Therefore, since a discontinuity in one or more 
[pjj+njj^ signals a IQPT, a discontinuity in Q (2,rt) also 
signals a IQPT. IfQ{2,n) is continuous and dxQ {2,n) 
shows a discontinuity or divergence, it signals a 2QPT. 
In this sense Q (2,n) is at least as good as the concur- 
rence/negativity to signal a QPT. Note that the previous 
result is valid only if the discontinuous/divergent quan- 
tities do not accidentally all vanish or cancel with other 
terms in Eqs. ^ and ^ (assumptions (b) and (c) in 



Ref. 01). An added bonus of our approach, however is 
that we do not need a further assumption, as in Ref. [7[ , 
related to the artificial/accidental divergences due to the 
maximization/minimization processes appearing in the 
definitions of the concurrence and the negativity. More- 
over, Q (2, n) is richer than the concurrence/negativity 
for signaling and classifying the order of a QPT, since it 
can be employed for the derivation of an EL, as we now 
demonstrate. 

We particularize our discussion to two-level (qubit) 
systems [ij. In this case G (2,n) is written as 



g(2,n) 



1 - 1 E 



a, 13=0 



(5) 



where erf, a = 1,2,3, are the Pauli operators and cr^ is 
the identity. Thus, as any measure dependent only on 
the two-particle reduced density matrix, G {2,n) is com- 
pletely determined by one- and two-point CFs. When- 
ever the system undergoes a second-order symmetry- 
breaking QPT, it will be reflected in one or more CFs 
and hence in the behavior of G{2,n). If the dominant 
(less rapidly decaying) CF decays with a power law 
((CT^^cr^") ~ n"'') at the critical point (implying a diverg- 
ing correlation length) and exponentially in its vicinity 



((cr^^CT^") ^ e "/'Sc;), so will Q{2,n) increase. For large 
n, g (2, n)7^g (2, oo) - (erf V^^"} V3. Hence, close to the 
critical point g (2, n) w ^ (2, oo) - Ce~2"/«c ^ where C is 
a constant, and g (2, n) increases exponentially fast, sat- 
urating for n We can then define an EL that is 
proportional to the correlation length, = The 
EL also diverges at the critical point with the same ex- 
ponent as ^c, such that ~ |A — Ad^"^. At A = Ac, 
for large n, g (2, n) ^ g (2, oo) — C'n^^''' , where C" is a 
constant, g {2,n) now increases as a power law with a 
power that is twice the CF exponent. Thus, g (2,n) in- 
herits all the universal properties of the CFs. Moreover, 
due to the g (2, n) scaling with n, at the critical point the 
entanglement is more distributed in the system (any two 
spins are entangled with the rest of the chain) than away 
from it, indicating ME Q prevails at the critical point. 
We emphasize that this result is quite general, applying 
to any collection of two-level systems with a second-order 
symmetry-breaking QPT. Next, we particularize to two 
specific cases in order to illustrate our general results. 

For an arbitrary spin-1/2 model presenting a IQPT, 
at least one of the CFs is discontinuous at the transition 
point. Thus, it is intuitive that g (2, n) is also discon- 
tinuous. A simple example is the frustrated two-leg spin- 
1/2 ladder discussed in Refs. 0, where all but the 
{afa"-), a = x,y,z, and (erf) expectation values vanish. 
The latter are discontinuous at the transition point but 
constant otherwise. The transition is clearly of first order 
and g (2, n) is able to signal it. 

The one- dimensional XY model in a transverse mag- 
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netic field is described by the following Hamiltonian 

^= - E{i [(1 + «+i + (1 - ^) ^F'^'+i] + '^^l ' 

(6) 

where N is the total number of spins (sites) and 7 > 
is the anisotropy. This Hamiltonian is symmetric under 
a global TT rotation about the z axis (c^^^' — > —c^'"^^), 
implying a zero magnetization in the x or y directions 
{{(Ji^^^) = 0). However, as the magnetic field h is de- 
creased (or J increased) this symmetry is spontaneously 
broken at A = J//i = 1 (in the thermodynamical limit) 
and a doubly-degenerate ground state with finite magne- 
tization (±i\f ) in the x direction develops, characterizing 
a ferromagnetic phase. It is possible then to define a sym- 
metric ground state (with (erf) = 0) as a superposition 
of these two degenerate ones. These states are of no use 
in practice, however, as they do not exist in real macro- 
scopic objects undergoing a phase transition ("clustering 
property"). We call non- symmetric or broken-symmetry 
states the ones in which there is a finite magnetization 
{{Sf = a'f 12) — ±A/). Note that at the paramagnetic 
phase there is no such distinction. By further decreasing 
the magnetic field, a second phase transition occurs at 
7^ -I- /i^ = 1. In this "third" phase the approach of the 
CFs to their saturation values is not monotonic but os- 
cillatory [3l ■ We should also say that this model reduces 
to the Ising model for 7 = 1, where only the first critical 
point exists, and to the XX model as 7 ^ 0. However, 
the XX model belongs to a different universality class 
and we consider here only < 7 < 1 [iG] . 

The XY model can be solved exactly and all the CFs 
are known [l^. To calculate t/(2,n) all we need is 
the one and two-point CFs (See Eq. (O). Due to 
the translational invariance of the model pij depends 
only on the distance n ~ \i — j\ between the spins and 
= (^i'^i+n) = Vn^- Remembering that is 
Hermitian and has a unitary trace we are left with nine 
independent elements of p^ , which may be functions of 
at most nine one- and two-point CFs. This number can 
be further reduced by the symmetries of the problem. 
The global symmetry under a tt rotation about the z 
axis yields {a<y^) = (afcrf ) = (erf erf) = (crfaf) = in 
the paramagnetic phase (A < 1). We end up with four 
elements: (o"^) and ((t"(t"), a = x,y,z. In the ferromag- 
netic phase (A > 1) this no longer holds since the Hamil- 
tonian symmetry is not preserved by the ground state 
and we need to evaluate the nine one and two-point CFs. 
The four CFs appearing in the paramagnetic phase and 
(af^^'') plus the three off-diagonal two-point ones were 
calculated in Ref. [ll. The first two pl'^ = pfy = for 
all values of 7 and A [16|. The last off-diagonal CF (p^^) 
was obtained exactly in terms of complex integrals whose 
calculation is cumbersome. However, we were able to ob- 
tain excellent bounds for it by imposing the positivity of 



the eigenvalues of pij 

With all the necessary CFs in hand Q (2, n) for 
the XY model reads 0(2, n) = l-i[2(CTp2 +2(ct|)2 + 

Fig. [1| we plot the lower and upper bounds for Q (2, 1) 
(by using the upper and lower bounds of p^^, respec- 
tively) as a function of A and for a few 7's. We first note 
that it is maximal at the critical point for any anisotropy 
(this is true throughout the interval < 7 < 1). Sec- 
ondly, the bounds obtained are very tight and can barely 
be distinguished for some anisotropics. Only in the fer- 
romagnetic phase and for 7 — > do the bounds become 
distinguishable. The derivative of (2, 1) with respect to 
A is depicted in Fig. [2] exhibiting, as expected, a diver- 
gence at the critical point. 
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Figure 1: (Color online) Upper and lower bounds of 5(2, 1) 
for the XY chain for three values of the anisotropy: 7 = 1 
(red/solid), 0.6 (blue/dashed), and 0.2 (black/dotted). 
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Figure 2: (Color online) Derivative of the lower bound of 
Q{2, 1) for three values of anisotropy: 7 = 1 (red/solid), 0.6 
(blue/dashed), 0.2 (black/dotted). The second phase transi- 
tion is also imprinted for the 7 = 0.2 as the curve crosses the 
abcissa at A = 1/ \/l — 7^. 

Now we analyse how Q{2, n) approaches its asymptotic 
value. It can be seen in Fig. [3]that Q (2, n) is an increas- 
ing function of the distance n. To study this behavior an- 
alytically we make use of the asymptotic form of the CFs 
of the XY model: for A < 1 [ig, (o-JcrJ+„) - n-^/^X^, 
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Figure 3: (Color online) Lower bound of Q (2, n) for 7 = 0.6 
and for three values of n: n — 1 (solid), 2 (dashed), and 7 
(long-dashed). 



) - n-3/2A«, and (cr|a|+„) - (a^)^ - n-2A2", 



with A2 = (1/A - VVA^ - (1 -7^))/(l - 7), while at 
the critical point (ctJctJ_^„) n^^/^, {a^^a^j^^) ^ jt'^I^, 



and (cr|cr|^_„) ~ (ct^)^ — We can see from these 

expressions that, for large values of n, the dominant cor- 
relation is, as expected, in the x direction. Thus, for 
large n we can write Q (2, ri) ^ Q (2, cxo) — (crJcrJ_^„)^/3, 
such that 



a (2, 71) ~ 5(2, 00) 
g(2,n) ~ 5(2, 00) 



A < Ac: 
A = Ac 



(7) 
(8) 



From these expressions, we see explicitly that, at the crit- 
ical point, the entanglement between two spins n sites 
apart increases as a power law of their distance, whereas 
away from the critical point it increases exponentially 
and saturates very fast. For the XY model the EL de- 
fined before reads (^e = 2{i-x) ' "^ti^re we have used that 
A2 ~ 1 + (A — l)/7 near the critical point. Note that 
diverges at the critical point as expected and that the 
ratio between and the correlation length is fixed: 
£,e/£,c = 1/2- Thus at the critical point the entangle- 
ment in the XY model is more distributed in the chain, 
as already indicated by the block entanglement Q . 

In conclusion, we related the non-analytic properties of 
the ground state energy to the non-analyticities of 5(2, n) 
for an arbitrary many-particle system. Thus, 5(2, n) is 
able to signal both the quantum phase transition (QPT) 
points and the order of the transtition. 5(2, n) is a mul- 
tipartite entanglement (ME) measure which, for many 
reasons (13i] . is operationally good. Since no maximiza- 
tion/minimization process is needed for its calculation, 
no accidental discontinuities or divergences will occur (in 
contrast to the concurrence or the negativity). Moreover, 
for two-level systems 5(2, rt) is simply related to one- 
and two-point correlation functions (CFs). Therefore, for 
those systems undergoing a second order QPT it is pos- 
sible to define a critical exponent and an entanglement 
length which is half the more familiar correlation length. 
We have exemplified those results with an explicit calcu- 



lation for the XY transverse field spin- 1/2 chain. This 
result adds strength to the conjecture by T. J. Osborne 
and M. A. Nielsen Q that at the critical point bipartite 
entanglement (as given by the concurrence/negativity) is 
not maximal due to entanglement sharing, since all the 
parties involved are entangled as the entanglement length 
diverges. In fact, what should be maximal and favored is 
the multipartite entanglement, as we have plenty demon- 
strated. It is worth mentioning that any knowledge of the 
behavior of ME can only be achieved via the generalized 
global entanglement (GGE) and not by any CF alone. 
We expect that these findings will contribute to the un- 
derstanding of the relevance of entanglement, specially 
ME, in QPTs. 

Note: After this work was completed we became aware 
of an independent derivation of the entanglement length 
for the XY model in terms of the two-site von Neumann 
entropy \17\. We point out, however, that the relatively 
simple form of 5(2, n), as given by the one and two-point 
CFs, allows it to be employed for the determination of 
the order of the QPT as well as for the derivation of 
an entanglement length for an arbitrary two-level system 
undergoing a second order QPT. 
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